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Toll-like receptors (TLRs) single nucleotide polymorphisms (SNPs) were analyzed in patients 
with papillary thyroid cancer (PTC; n = 133) and their clinicopathologic features and age- 
matched controls (n = 321) using direct sequencing. PTC patients were divided into 
subgroups according to size, number, location, extrathyroidal invasion and lymph node 
metastasis. The two SNPs of TLR2 gene were not associated with the development of PTC. 
In clinical analysis, two SNPs were associated with location of cancer (rs3804099, 
P= 0.032, OR, 0.52; 95% CI, 0.28-0.96 in log-additive model; rs3804100, P= 0.039, 
OR, 0.46, 95% CI, 0.22-0.96 in codominanti model; P= 0.018, OR, 0.42, 95% CI, 0.21- 
0.87 in dominant model; P= 0.011, OR, 0.46, 95% CI, 0.25-0.85 in log-additive model). 
The allele frequencies of two SNPs also showed significant associations with location of 
cancer (rs3804O99, P= 0.046, OR, 0.57, 95% CI, 0.33-0.99 and rs3804100, P= 0.019, 
OR = 0.52, 95% CI = 0.30-0.90). However, two SNPs were not associated with the 
clinicopathologic features of PTC. It is suggested that 71^2 polymorphisms may contribute 
to the clinicopathologic features of PTC, especially the PTC in both lobes. 

Keywords: Papillary Thyroid Cancer; Single Nucleotide Polymorphisms; Clinicopathologic 
Characteristics; Toll-Like Receptors 
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INTRODUTION 

Thyroid cancer is characterized by the most common endocrine 
malignant carcinoma and is one of the fastest growing cancer 
diagnoses worldwide. There are four major types of thyroid can- 
cers that are clinically significant: papiUary, follicular, anaplas- 
tic, and medullary thyroid cancers. Of these, papillary thyroid 
cancer (PTC) is the most common type, showed about 80% of 
all cases of thyroid cancer It is well-known that ionizing radia- 
tion exposure is a main risk factor (1). However, the exact etiol- 
ogy of PTC has been still unknown. 

The Toll-like receptors (TLRs) family of transmembrane re- 
ceptors plays a fundamental role in pathogen recognition and 
activation of innate immunity. Ligation of TLRs results in sub- 
stantial immune responses that may be directed against tumor- 
associated antigens (2). A latent role for radiation in signaling 
'danger' and, perhaps, in the activation of antigen presenting 



cells on immune system has been suggested (3-6). Radiation 
may also activate effectors of innate immunity through TLR- 
dependent mechanisms, thereby up-regulate the adaptive im- 
mune response to cancer (2). The TLR2 gene is located in the 
chromosome 4p32 protein coding gene, and is expressed most 
abundantly in peripheral blood leukocytes, and mediates host 
response (http://www.ncbi.nlm.nih.gov/gene). 

Hashimoto's thyroiditis is the most common gradual auto- 
immune-mediated thyroid failure and it is developed goiter at 
times. Some studies have found an association between Hashi- 
moto's thyroiditis and PTC (7). Previous study was demonstrat- 
ed that Hashimoto's thyroiditis was associated with an increased 
risk of developing PTC (7). Also, genetic factors seem to be in- 
volved as a risk factor for PTC (8-12). Indeed, several single nu- 
cleotide polymorphisms (SNPs) are thought to contribute to 
susceptibility to PTC (13-15). 

In this study, we investigated whether synonymous SNPs in 



© 2012 The Korean Academy of Medical Sciences. 

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.0rg/licenses/by-nc/3.o) 
which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited. 



pISSN 1011-8934 
elSSN 1598-6357 



JKMS 



Kim MK, etal. • Association ofTLR2 Polymorphism with Papillary Thyroid Cancer 



TLR2 contribute to the development of PTC. We also assessed 
the relationships between TLR2 SNPs and the cUnicopathologic 
characteristics of PTC. 

MATERIALS AND METHODS 
Patients and controls 

Study subjects consisted of PTC patients (n = 133, 48 males and 
85 females) and controls (n = 321, 129 males and 192 females). 
The mean ages of the PTC and control groups were 54.7 ± 12.3 
yr (mean ± SD) and 56.3 ± 11.9, respectively (Table 1). PTC pa- 
tients were recruited among participants visiting Kyung Hee 
University Medical Center, Seoul, Republic of Korea. Control 
subjects were enrolled from healthy participants examined in a 
general health check-up program. Participants with cancers, 
thyroid diseases, or any other severe diseases were excluded. 
PTC diagnosis was confirmed by pathologic examination. Pa- 
tients with anaplastic carcinoma, follicular carcinoma, double 
primary of PTC and follicular carcinoma, follicular variant of 
PTC, or nodular hyperplasia were excluded. 

Patient subgroups 

To assess the relationship between TLR2 SNPs and the clinical 
pathologic characteristics of PTC, patients were divided into 
subgroups according to the size (< 1 cm and > 1 cm), number 
(unifocality and multifocality), location (one lobe and both 



lobes), extrathyroidal invasion (present and absent), and lymph 
node metastasis (present and absent). Demographic features of 
PTC patients are summarized in Table 1. 

SNP selection and genotyping 

For the selection among TLR2 SNPs, we searched the synony- 
mous SNPs of the TLR2 gene in the SNP database of the National 
Center for Biotechnology Information (http://www.ncbi.nlm. 
nih.gov/SNP BUILD132). Because 16 missense SNPs of the TLR2 
gene had an unknown heterozygosity or minor allele frequency 
(MAP) below 0.05, all missense SNPs were excluded. Out of 13 
synonymous SNPs of the TLR2 gene, there were 1 1 SNPs with 
MAF below 0.05. Finally, two SNPs (rs3804099, Asnl99Asn; 
rs3804100, Serl32Ser) were selected. Blood samples for DNA 
extraction from all subjects were collected. Genomic DNA was 
extracted with 200 |iL whole blood using DNA Isolation Kit for 
Cells and Tissues (Roche, Indianapolis, IN, USA) and stored at 
-20°C before use. SNP genotyping was determined by direct se- 
quencing. Polymerase chain reactions (PCRs) were performed 
using the primers for two synonymous SNPs (Table 2). PCR com- 
prised 40 cycles at 94°C for 30 sec, 58°C for 30 sec, 72°C for 30 sec, 
and 1 cycle at 72°C for 5 min for the final reaction. The PCR prod- 
ucts were sequenced by an ABI PRISM 3730XL analyzer (PE 
Applied Biosystems, Foster City, CA, USA). Sequencing data 
were analyzed using SeqManll software (DNASTAR, Madison, 
WI, USA). 



Table 1. Clinical cinaracteristics of the study population 



Parameters 


PTC (No. = 133) 


Control (No. = 321) 


Age (mean + SD) 


55.7 + 12.8 


56.3 ±11.9 


Sex (male/female, No.) 


48/85 


129/192 


Size of cancer (No.) 
< 1 cm 
> 1 cm 


61 

70 




Number of cancer (No.) 
One 
IVIultiple 


83 
45 




Location of cancer (No.) 
One lobe 
Both lobes 


63 
65 




Extrathyrodial invasion (No.) 
Absent 
Present 


67 
62 




Cervical lymph node metastasis (No.) 
Absent 
Present 


84 
40 





PTC patients having inappropriate clinical data were excluded. Because of PTC pa- 
tients having inappropriate clinical data, the total number of PTC phenotype are dif- 
ferent. PTC, papillary thyroid cancer; No., number of subjects; SD, standard deviation. 



Statistical analysis 

Hardy- Weinberg equilibrium (HWE) was estimated using SNP- 
Stats (http://bioinfo.iconcologia.net/index.php7module = Snp- 
stats) in both the patient and control groups. Multiple logistic 
regression models were applied to obtain odds ratios (ORs), 95% 
confidence intervals (CIs), and P values (codominantl, codom- 
inant2, dominant, recessive, and log-additive models). When 
the numbers of genotype and allele were less than 5, Fisher's 
exact test was performed. Data analysis was performed using 
SPSS 18.0 (SPSS, Chicago, IL, USA) and SNPstats (http://bioin- 
fo.iconcologia.net/index.php?module = Snpstats). Linkage dis- 
equilibrium (LD) block and haplotypes were evaluated using 
Haploview version 4.2 (Daly Lab Inc., Cambridge, MA, USA). 
P < 0.05 was considered significant. 

Ethics statement 

This study was approved by the institutional review board of the 
Medical Research Institute, Kyung Hee University Medical Cen- 



Table 2. Primer sequences for each SNP 



SNP 


Forward (5 '-3') 


Reverse (5 '-3') 


Size (bp) 


rs3804099 


ATTGCAAATCCTGAGAGTGGGA 


GTTTCACCAGTGGATAGTTCTG 


291 


rs3804100 


GAAAAATTCAGCCTGTGAGGAT 


GGAGGCATCTGGTAGAGTCATC 


358 



SNP single nucleotide polymorphism; bp, base pair 
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ter, Seoul, Korea (20040915). Written informed consent was ob- 
tained from all subjects. 

RESULTS 

We genotyped two synonymous SNPs in the TLR2 gene. The 
genotype and allele frequencies of two synonymous SNPs are 
presented in Table 3. Multiple logistic regression analysis with 
adjustment for age and gender was performed: codominantl 
(major allele homozygotes vs heterozygotes), codominant2 (ma- 
jor allele homozygotes vs minor allele homozygotes), dominant 
(major allele homozygotes vs heterozygotes + minor allele ho- 
mozygotes), recessive (major allele homozygotes + heterozy- 
gotes vs minor allele homozygotes), and log- additive (major al- 
lele homozygotes vs heterozygotes vs minor allele homozygotes). 



The two SNPs were in Hardy- Weinberg equilibrium in the con- 
trol group (rs3804099, P = 0.30; rs3804100, P = 0.11). Our two 
synonymous SNP in Table 3 (rs3804099, Asnl99Asn; rs3804100, 
Serl32Ser) of TLR2 was not associated between control and 
PTC patients (rs3804099, P = 0.680 in allele [reference T vs C]; 
rs3804100, P = 0.850 in allele [reference Tvs C]). 

Next, we assessed the relationship between TLR2 SNPs and 
the clinical pathologic features of PTC. In the location of cancer 
(one lobe vs both lobes), the SNP (rs3804099) was related to lo- 
cation of PTC in genotype (rs3804099, P = 0.032, OR, 0.52; 95% 
CI ,0.28-0.96 in log-additive model [reference T/T vs T/C vs C/ 
C]) and allele distributions (P = 0.046, OR, 0.57; 95% CI = 0.33- 
0.99, reference T vs C). Another SNP (rs3804100) was signifi- 
cantly associated with location of PTC in genotype (rs3804100, 
P = 0.039, OR, 0.46, 95% CI, 0.22-0.96 in codominantl model 



Table 3. Frequencies of genotype and allele of TLR2m control and patients with papillary thyroid cancer (PTC) 



SNP 




Type 


Control 


PTC 


Model 


OR (95% CI) 


P value 

-al 




n (%) 


n (%) 


rs3804099 


Genotype 


T/T 


157 (48.9) 


61 (46.9) 


Codominantl 


1.23 (0.80-1.88) 


0.350 


Asn199Asn 




T/C 


129 (40.2) 


61 (46.9) 


Codominant2 


0.56 (0.24-1.28) 


0.170 






C/C 


35 (10.9) 


8 (6.2) 


Dominant 
Recessive 
Log-additive 


1.08 (0.72-1.62) 
0.51 (0.23-1.13) 
0.93 (0.68-1 .27) 


0.710 
0.080 
0.630 




Allele 


T 

C 


443 (69.0) 
199(31.0) 


183(70.4) 
77 (29.6) 




1 

0.94 (0.68-1.28) 


0.680 


rs3804100 


Genotype 


T/T 


165(51.4) 


64(48.1) 


Codominantl 


1.30 (0.85-1.99) 


0.220 


Ser450Ser 




T/C 


122 (38.0) 


61 (45.9) 


Codominant2 


0.58 (0.25-1.33) 


0.200 






C/C 


34 (10.6) 


8 (6.0) 


Dominant 
Recessive 
Log-additive 


1.14(0.76-1.71) 

0.089 
0.96 (0.71-1.31) 


0.530 
0.810 




Allele 


T 


452 (70.4) 


189(71.1) 




1 








C 


190 (29.6) 


77 (28.9) 




0.97(0.71-1 .33) 


0.850 


The P values were calculated from logistic regression analyses adjusting sex and age (codominantl , A/A vs A/B; codominant2, A/A vs B/B; dominant A/A vs A/B+B/B; reces- 


sive, A/A+A/B vs B/B; log-additive, A/A vs A/B vs B/[ 


3). Missing data of genotype were excluded for exact analysis. PTC, papillary thyroid cancer; SNP, single nucleotide poly- 


morphism; OR, odds ratio; CI, confidence interval. 












Table 4. Genotype and allele frequencies of SNPs of 7I/?2gene in PTC patients with one lobe and PTC patients with both lobes 






SNP 




Type 


One lobe 


Both lobes 


Model 


OR (95% CI) 


P value ^ 




n (%) 


n (%) 


rs3804099 


Genotype 


T/T 


23 (37.1) 


34 (54.0) 


Codominantl 


0.53 (0.25-1.13) 


0.100 


Asn199Asn 




T/C 


33 (53.2) 


27 (42.9) 


Codominant2 


0.24 (0.04-1.34) 


0.100 






C/C 


6 (9.7) 


2 (3.2) 


Dominant 
Recessive 
Log-additive 


0.49 (0.24-1.02) 
0.34 (0.06-1.78) 
0.52 (0.28-0.96) 


0.054 
0.170 
0.032 




Allele 


T 


79 (63.7) 


95 (75.4) 




1 








C 


45 (36.3) 


31 (24.6) 




0.57 (0.33-0.99) 


0.046 


rs3804100 


Genotype 


T/T 


23 (36.5) 


37 (56.9) 


Codominantl 


0.46 (0.22-0.96) 


0.039 


Ser450Ser 




T/C 


34 (54.0) 


26 (40.0) 


Codominant2 


0.21 (0.04-1.18) 


0.080 






C/C 


6 (9.5) 


2(3.1) 


Dominant 
Recessive 
Log-additive 


0.42 (0.21-0.87) 
0.32 (0.06-1 .68) 
0.46 (0.25-0.85) 


0.018 

0.150 
0.011 




Allele 


T 


80 (63.5) 


100 (76.9) 




1 








C 


46 (36.5) 


30 (23.1) 




0.52 (0.30-0.90) 


0.019 



The P values were calculated from logistic regression analyses adjusting sex and age (codominantl , A/A vs A/B; codominant2, A/A vs B/B; dominant A/A vs /VB-rB/B; reces- 
sive, A/A-rA/B vs B/B; log-additive, A/A vs A/B vs B/B). Bold numbers mean significance association. Missing data of genotype were excluded for exact analysis. PTC, papillary 
thyroid cancer; SNR single nucleotide polymorphism; OR, odds ratio; CI, confidence interval. 
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Table 5. Genotype and allele frequencies of SNPs of TLR2qene in PTC patients with tumor size < 1 cm and PTC patients with tumor size > 1 cm 









< 1 cm 


> 1 cm 








SNP 




Type 


(Tumor size) 


(Tumor size) 


Model 


OR (95% CI) 


P value \ 








n (%) 


n (%) 








rsooU4uyy 


Genotype 


T/T 
1/ 1 




OO (0 1 .0) 


Codominanti 


n c;o /n OQ i or;\ 

u.oy (u.zo- 1 .^loj 


n i 7n 
U. 1 /U 


Msn 1 yyAsn 




T/r 

1 /Lj 


o 1 p 1 . / ) 


on lAO R\ 


Codominant2 


U.DO (U. 1 o-o. 1 o) 


n KQn 






U/U 


/I la 7\ 
4 (b./) 


4 (o.y) 


Dominant 


n cn /n on i '^a\ 

u.du (u.zy- 1 .£14) 


n i 7n 
U. 1 /U 












Recessive 


n on i{\ on Q GA\ 

U.oy (u.zu-o.y4j 


n RHn 












Log-additive 


n fiQ /n ofi 1 oR\ 
u.oy (u.oo- 1 .ZDj 


n OQn 




Allele 


T 


81 (67.5) 


99 (72.8) 




1 








C 


39 (32.5) 


37 (27.2) 




0.78 (0.45-1 .33) 


0.360 


rs3804100 


Genotype 


T/T 


27 (44.3) 


36(51.4) 


Codominanti 


0.62 (0.30-1 .30) 


0.200 


Ser450Ser 




T/C 


31 (50.8) 


29(41.4) 


Codominant2 


1.18(0.25-5.59) 


0.840 






G/C 


3 (4.9) 


5 (7.1) 


Dominant 


0.67 (0.33-1.36) 


0.270 












Recessive 


1 .49 (0.33-6.74) 


0.600 












Log-additive 


0.81 (0.45-1 .45) 


0.480 




Allele 


T 


85 (69.7) 


101 (72.1) 




1 








C 


37 (30.3) 


39 (27.9) 




0.89 (0.52-1.51) 


0.660 


The P vaiues were caiculated from iogistic regression analyses adjusting sex and age (codominanti , A/A vs A/B; codominant2, A/A vs B/B; dominant /VA vs A/B+B/B; reces- 


sive, A/A+ZVB vs B/B; log-additive, A/A vs A/B vs B/B). Missing data of genotype were excluded for exact analysis. PTC, papillary thyroid cancer; SNP single nucleotide poly- 


morphism; OR, odds ratio; CI, confidence interval. 












Table 6. Genotype and allele frequencies of SNPs of 7I/?2gene in PTC patients with unifocality and PTC patients with multifocality 












Unifocality 


Multifocality 








SNP 




Type 


Model 


OR (95% CI) 


P value ^ 






n (%) 


n (%) 






rs3804099 


Genotype 


T/T 


34 (42.5) 


23(51.1) 


Codominanti 


0.75 (0.35-1.61) 


0.470 


Asn199Asn 




T/C 


40 (50.0) 


20 (44.4) 


Codominant2 


0.45 (0.08-2.49) 


0.360 






C/C 


6 (7.5) 


2 (4.4) 


Dominant 


0.71 (0.34-1.50) 


0.370 












Recessive 


0.51 (0.10-2.76) 


0.420 












Log-additive 


0.71 (0.38-1.33) 


0.290 




Allele 


T 


108 (67.5) 


66 (73.3) 




1 








C 


52 (32.5) 


24 (26.7) 




0.76 (0.43-1 .34) 


0.340 


rs3804100 


Genotype 


T/T 


36 (43.4) 


24 (53.3) 


Codominanti 


0.71 (0.33-1.52) 


0.380 


Ser450Ser 




T/C 


41 (49.4) 


19(42.2) 


Codominant2 


0.47 (0.09-2.60) 


0.400 






C/C 


6 (7.2) 


2 (4.4) 


Dominant 


0.68 (0.33-1 .42) 


0.310 












Recessive 


0.56 (0.11-2.96) 


0.480 












Log-additive 


0.70 (0.38-1.31) 


0.260 




Allele 


T 


113(68.1) 


67 (74.4) 




1 








C 


53(31.9) 


23 (25.6) 




0.73 (0.41-1.30) 


0.290 



The P values were calculated from logistic regression analyses adjusting sex and age (codominanti , /VA vs A/B; codominant2, A/A vs B/B; dominant /VA vs A/B-rB/B; reces- 
sive, A/A-rA/B vs B/B; log-additive, A/A vs A/B vs B/B). Missing data of genotype were excluded for exact analysis. PTC, papillary thyroid cancer; SNP single nucleotide poly- 
morphism; OR, odds ratio; CI, confidence interval. 



[reference T/T vs T/C]; P = 0.018, OR, 0.42, 95% CI, 0.21-0.87 in 
dominant model [reference T/T vs T/C -i- C/C]; P = 0.011, OR, 
0.46; 95% CI, 0.25-0.85 in log-additive model [reference T/T vs 
T/C vs C/C], respectively) and allele distributions (P - 0.019, 
OR, 0.52; 95% CI, 0.30-0.90, reference T vs C) (Table 4). Howev- 
er, no association was evident with clinical pathologic features 
ofPTC (Tables 5, 6). 

Linkage disequilibrium (LD) block between rs3804099 and 
rs3804100 was determined using Haploview version 4.2. Since 
the LD block was made (D' = 1.00 andr^ = 0.94, data not shown), 
the analysis of haplotypes was performed. The frequency of TT 
was 0.69, CC was 0.297 and CT was 0.013, however, the associa- 
tion was not observed (data not shown). We calculated sample 
power to verify our data (http://pngu.mgh.harvard.edu/~pur- 
ceU/gpc/cc2.htnil). In this study, the sample power of each SNP 

1336 http;//jl<ms.org 



was 0.745 (rs3804099, alpha level = 0.05, the number of cases for 
70% power = 91) and 0.730 (rs3804100, alpha level = 0.05, the 
number of cases for 70% power = 95). These sample powers were 
sufficient to detect positive associations between SNPs and PTC. 

DISCUSSION 

The functions of TLRs have provided potential insight into the 
cancer development. TLR-dependent mechanisms contribute 
to activate effectors of innate, thereby augmenting the adaptive 
immune response to cancer (2, 16, 17). 

Numerous studies have investigated TLR2 polymorphisms 
in different populations and different diseases with conflicting 
results. The two presently-studied TLR2 SNPs, rs3804099 and 
rs3804100, were reportedly associated with sarcoidosis with 
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cutaneous manifestations in a Japanese population (18), and 
rs3804099 was suggested as a relevant risk estimate for the de- 
velopment of sepsis and multiple organ dysfunction in Chinese 
patients with major trauma (19). In addition, TLR2 -196 to -173 
del, +597 T > C and +1350 T > C polymorphisms have been as- 
sociated with asymptomatic bancroftian fQariasis in Thailand 
(20). In a TLR2-related cancer study, SNPs of TLR2 were not as- 
sociated with colorectal cancer in individuals of three races (20 
Malays, 20 Chinese and 18 Indians) (21). On the contrary, the 
association with microsateUite GT polymorphisms of TLR2 gene 
and sporadic colorectal cancer among Croatians has been re- 
ported (22). 

Here, we sought to determine the relationships between the 
TLR2 SNPs and PTC in a Korean population. We also assessed 
the relationships between TLR2 SNPs and the cUnicopathologic 
characteristics of PTC. The findings of our study suggest that no 
significant differences exist in the frequency of TLR2 genotypes 
and alleles in the PTC cases comparison with the controls. How- 
ever, TLR2 was associated with clinicopathologic features, es- 
pecially in tiie locations of PTC. The T/T frequency of rs3804099 
was different between the one lobe and both lobe groups (37.1% 
vs 54.0%, respectively). The T allele frequencies of rs3804099 in 
PTC patients with involvement of both lobes (75.4%) were high- 
er than those in PTC patients with one lobe involvement (63.7%). 
The T allele frequencies of rs3804100 in PTC patients witii in- 
volvement of both lobes (76.9%) were higher than those in PTC 
patients with involvement of one lobe (63.5%). Therefore, the T 
allele of rs3804099 and rs3804100 may be correlated with the 
location of PTC in the Korean population. Thyroid lobectomy 
alone may be sufficient tteatment for smaU ( < 1 cm), low- risk, 
unifocal, intrathyroidal papillary carcinomas in the absence of 
prior head and neck irradiation or radiologicaUy or clinically in- 
volved cervical nodal metastases. But patients who have bilat- 
eral PTC should undergo total thyroidectomy. Therefore, tumor 
bilaterality can be used to determine surgical extent (23). The 
relationship between TLR2 polymorphisms and cancer loca- 
tion is significant. 

This study has some limitations. The sample size of patients 
is small and control subjects did not receive the thyroid ultra 
sonography. Considering that the incidence of thyroid cancer 
(0.5-10/100,000 persons) is relatively low, more studies with 
larger numbers of patients are needed to verify our results. 

In conclusion, based on our case-control association study of 
SNPs in TLR2 genes in patients with PTC and conttol subjects, 
significant associations are reported between polymorphisms 
of the TLR2 gene and PTC in both lobes. The results suggest that 
the TLR2 polymorphisms may be associated with the clinico- 
pathologic features of PTC in the Korean population, especially 
concerning the location of cancer. 
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